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1. Two dimensional SU(N) X SU(N) principal chiral model 

The two dimensional SU(N) X SU(N) principal chiral model is similar to four dimensional 
SU(N) gauge theory in many respects [jl|]. The continuum action is given by 

S=j f d 2 xTrd^g(x)d^(x) (1.1) 

where g(x) E SU(N). The global symmetry group SU(N)^x SU(N)# reduces down to a single 
SU(N) "diagonal subgroup" if we make a translation breaking "gauge choice", g(0) = 1. This 
model is asymptotically free and there are N — 1 particle states with masses 

sinf— ) 

M R = M— ±*LL, l<R<N-\. (1.2) 
sm (f) 

The states corresponding to the R-th mass are a multiplet transforming as an R component antisym- 
metric tensor of the diagonal symmetry group. 

The two point function W = g(0)g^(x) plays the role of Wilson loop with the separation x 
playing the role of area. We expect the behavior to be perturbative for small x. On the other hand, 
non-perturbative effects become important for large x. 

One expects 

Gr(x) = W«(0)/W))~QQ e -*H (1.3) 

where Xr is the trace in the /^-antisymmetric representation. Comparison with the heat-kernel 
representation of the characteristic polynomial associated with the Wilson loop operator in two 
dimensional large N QCD [Q] suggests the following connections: 

• The two point correlator, W(d) = g(0)g^(d), is analogous to the Wilson loop operator. 

• M\x\ is analogous to the dimensionless area, t. 

Based on this analogy, we hypothesize ^ that the characteristic polynomial, det (z — g(0)g'(d)), 
will undergo a transition at some value d c . The universal behavior at this transition will be in the 
same universality class as two dimensional large ,/V QCD. 

2. Setting the scale 

Numerical measurement of the correlation length using the lattice action 

S L = -2Nb^Tr[g(x)g^x + ^)} (2.1) 



and 

2 = 1 I,y* 2 Gi(x) 

q ° 4 £,Gi(*) 
yields the following continuum result [Q] : 



(2.2) 



M& = 0.991(1) (2.3) 
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We use to set the scale and it is well described by 

£ G = 0.991 



2-n 

e 4 



167T 



VE exp (|) (2.4) 



in the range 1 1 < £ G < 20 with 



r /„x t,M ( 1 1 0.000545 0.00095 0.00043 

E = 1 - ^(TVfe^'O)]} = ^ + 5^ + - j, + -p- 

The above equations will be used to find a b for a given £ . 

3. Smeared SU(N) matrices 

Well defined operators are obtained using smeared matrices. We start with g(x) = go(x) and 
one smearing step takes us from g t (x) to g t+ \ (x) using the following procedure. Define Z t+ \ (jc) by: 

Z t+l (x)=^[gHx)g t (x + ^-\]. (3.1) 
±n 

Construct anti-hermitian traceless SU (N) matrices A t+ \ (x) 

A t+1 (x) = Z t+l (x) - Zj +l (x) - ^Tr(Z t+1 (x) - z\ +l (x)) = -A t t +1 (x). (3.2) 

Set 

L t+ i(x) =exp[/A ?+ i(x)]. (3.3) 
g t+ \ (x) is defined in terms of L t+ \ (x) by: 

8t+i(x) = g t (x)L t+ i(x). (3.4) 

This procedure is iterated till we reach g n {x) and the smearing parameter is defined by % = nf. 
For a fixed E,q, the parameter t is fixed such that t/<^ g remains unchanged. We set « = 30 in our 
numerical simulations and this was found sufficiently large to eliminate a dependence on the two 
factors, / and n, individually. 

4. Numerical details 

We need L/E,g>1 to minimize finite volume effects. We worked in the range 11 < E, G < 20 
and therefore we chose L = 150. We used a combination of Metropolis and over-relaxation at each 
site x for our updates. The full SU(N) group was explored. 200-250 passes of the whole lattices 
were sufficient to thermalize starting from g(x) = 1. 50 passes per step were enough to equilibrate 
if <§g was increased in steps of 1. 

The test of the universality hypothesis proceeds in the same manner as for the three dimen- 
sional large N gauge theory. We defined the characteristic polynomial, F(y,d), as 

F(y,d) = (det(e y/2 + e- y/2 W(d))} (4.1) 
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Figure 1: Behavior of Q. as a function of a in the scaling region. 



We perform a Taylor expansion, 

F(y,d,N)= C (d,N) + C 2 {d,N)y 2 +C 4 (d,N)y 4 + ... 
since F(y,d) is an even function of y. It is useful to define 

C (d,N)C 4 (d,N) 



a(d,N) 



C%(d,N) 



(4.2) 



(4.3) 



1 for short distances d < d c and £1 = ^ for 
oo in the vicinity of d = d c using 



which resembles a Binder cumulant. 

As N — > oo, n(<i,oo) is a step function with Q. - ^ 
long distances, d > d c . Zooming in on the step function as N 
the scaling variable a c< \/N(d — d c ), we obtain Fig. [jj. 

We use £l(a = 0) = 0.364739936 to obtain the critical size d c in the following manner. Given 
an N and a we find the d c that makes the Binder cumulant Q.(d c ,N) = 0.364739936 as shown 
in Fig. ||. We look at d c as a function of § for a given A^. This gives us the continuum value of 
d c /% for that Af. This extrapolation is shown in Fig. || for N = 30. We then take the large Af limit 
as shown in Fig. ^| and it gives us 

^|iv=oo = 0.885(3) (4.4) 



Further substantiation of the universal behavior can be given by comparing the eigenvalues 
distribution in the model to the Durhuus-Olesen eigenvalue distributions in two dimensional QCD. 
This is shown for one example each on either side of the critical point in Fig. ^ and very close to 
the critical point in Fig. 0. We use 2k = t to match with the notation in [S]. 



4 



Large N transition in the 2D SU(N)xSU(N) nonlinear sigma model. 



Rajamani Narayanan 



0.05- 



a 



-0.05 



-0.1 



-0.15 





■ ••7 ♦ 

T- ~'„'-H.,,, 






N=15 spline 


o 


N=15data 




- N=20 spline 


A 


N=20 data 




-N=25 spline 


V 


N=25 data 




N=30 spline 


□ 


N=30 data 




— N=40 spline 




N=40 data 



-0.2 
0.45 



0.5 



0.55 



0.6 



0.65 



0.7 



0.75 



0.8 



0.85 



Figure 2: Plot of Q(d) after the subtraction of Q.(a = 0) = 0.364739936 as a function of d/| G . 
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Figure 3: Extrapolation to continuum of d c j for N = 30. 
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Figure 4: Extrapolation of the continuum d c j % to infinite N. 
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Figure 6: An example of an almost critical eigenvalue distribution. 
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